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Abstract. The recent measurements of the electric form factor of the neutron suggest that its shape may be interpreted 
as a smooth broad distribution with a bump at Q 2 ~ 0.3 (GeV/c) 2 superimposed. As a consequence the corresponding 
charge distribution in the Breit frame shows a negative charge extending as far out as 2 fm. It is natural to identify this 
charge with the pion cloud. This realisation is then used to reanalyse all old and new data of the electric and magnetic 
from factors of the proton and the neutron by a phenomenological fit and by a fit based on the constituent quark model. 
It is shown that it is possible to fit all form factors coherently with both ansaetzen and that they all show the signal of 
the pion cloud. 

PACS. 14.20.Dh - 13.40.Gp - 21.10.Ft 



1 Introduction 

Form factors encode unique information about the internal struc- 
ture of a scatterer provided they are determined with sufficient 
precision over a sufficiently large range of momentum transfer. 
Depending on the interaction, the Fourier transformation of the 
form factors gives the spatial distribution of e.g. mass, charge 
or magnetisation, which provides insight into several aspects of 
the internal structure of the scatterer: 

- the constituents present in the system, 

- their interaction, 

- and their wave functions. 

Therefore, form factors represent very significant tests of any 
model of the scatterer. 

The nucleon is realized in nature in two species, the pro- 
ton with one charge unit and the neutron with no net charge. 
While the proton should dominantly be describable by the s- 
state wave functions of the two up and the one down constituent 
quark, these contributions cancel to first order in the neutron 
and its electric form factor should be zero in this approxima- 
tion. In this simple picture it is assumed that the quarks are 
dressed by gluons and sea quarks forming "constituent quarks" 
which represent effective Fermions with equal masses of about 
one third of the nucleon mass. However, already before the re- 
alization of the quark-gluon structure of the nucleon the per- 
ception of a pion cloud around the nucleon was used in order to 
account for the Yukawa interaction between the nucleons. After 
early pure quark models of the nucleon, like the MIT bag, the 
necessity of a pion cloud was soon realized in order to preserve 
chiral symmetry at the nucleon surface, and the "little bag" and 
"cloudy bag" models were invented. Based on the fundamen- 
tal chiral symmetry of the QCD Lagrangian "chiral dynamics" 
was developed identifying the pion as the almost Goldstone 



Boson of the strong interaction. The "chiral perturbation the- 
ory" based on it has shown through many experimental tests 
that indeed the pion is a decisive constituent of the nucleon be- 
sides the elementary quarks and gluons. It is the purpose of this 
paper to check whether and how the pion cloud is reflected in 
the nucleon form factors. 

Because of its zero charge, the contribution of the bare neu- 
tron n° to the electric form factor of the physical neutron n 
is small, and the dissociation of a nucleon into its counterpart 
(here: the proton), and a charged pion (here: a negative pion) 
should emerge most clearly in the neutron's electric form fac- 
tor Geu- We therefore start out in section 2 with a discussion 
of Geu for which now data exist from polarisation measure- 
ments [1-14] having a smaller model dependence than previ- 
ous determinations. In section 3 we give an overview of the 
existing relevant form factor measurements and of our data se- 
lection for the present investigation. In section 4 we show that, 
at the percentage level, the peculiar structure observed in Gb„, 
namely a kind of bump around Q 2 w 0.2 — 0.3 (GeV/c) 2 , is 
also present in the other form factors Ge p , Gm p , and Gmu- 
While this is discussed in section 4 in terms of a purely phe- 
nomenological ansatz for the form factors, we show in section 5 
how this can be viewed in the light of the decomposition of the 
nucleon states into a constituent quark core and a polarisation 
term reflecting the contribution of the pion cloud. The findings 
are discussed in the concluding section 6. 



2 The triggering conjecture: The electric 
form factor of the neutron G En 

Previous efforts to determine the electric form factor of the 
neutron, e. g. from elastic electron scattering on the deuteron, 
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were hampered by severe model-dependencies of the results, 
which therefore were uncertain to about 50%. The emerging 
results were describable by the so called Galster parameterisa- 
tion, which started out from the usual dipole fit, which repro- 
duced Ge p , Gm p , and Gmu reasonably well and which was 
multiplied by some appropriate function in order to account for 
the condition GEn{Q 2 = 0) = required by the vanishing 
charge of the neutron. This Galster form is given by 



G E n{Q 2 ) = 



a G t 



1 



(1 + b G t) 



(l + Q 2 /m 2 D ) 2 



(1) 



where r = Q 2 /(2m n ) 2 and m n = 0.939 GeV/c 2 is the neu- 
tron mass. The parameter m 2 D was taken as the standard dipole 
value m 2 D = 0.71 (GeV/c) 2 and ag = 1.73 in order to repro- 
duce the measured root mean square radius of the neutron of 
(r 2 ) = -6 dG E n(Q 2 )/dQ 2 \ Q 2 =0 = -0.115 fm 2 as deter- 
mined from the scattering of thermal neutrons [15]. Thus the 
only parameter free to be fitted to the data was bo, and it was 
determined to b<j = 4.59. The Galster form has no particu- 
lar theoretical justification and may rather hide the essential 
physics. 

The collected data for Geu determined recently from po- 
larisation measurements are depicted in fig. 1. These 15 data 
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Fig. 1. The Geu data from polarisation measurements. The coding 
for reactions with the deuteron as a neutron target is: open square [1], 
filled diamonds [2,3], open diamond [4], open star [9], open triangle 
[14], open pentagon [12, 13], and filled triangle [10, 1 1], the measure- 
ments with 3 He are shown as filled squares [5-8]. The full curve de- 
picts the fit of the parameters of eq. (2) to the data, the dashed dotted 
curve is a variant with slightly changed parameters as explained in the 
text, while the dotted curve is a fit using the Galster form, i.e. eq. (1). 



points, which are not hampered by model assumptions, have 
been taken with 8 very different experimental setups, and the 
data points taken with the same setups were taken over periods 
separated by long time intervals. Also, the setups had very dif- 
ferent systematic errors and corrections due to nuclear binding 
effects. Therefore, it is justified to consider the data as statisti- 
cally independent. Since the corrections are less certain for the 
measurements on 3 He than for the loosely bound deuterium, 
the measurements on the two targets are distinguished in fig. 1 
by markedly different symbols. It is not the aim of this paper, 



however, to discuss critically these experiments but just to take 
this data set seriously and to investigate its essential features. 

It is evident from fig. 1 that the data can be as well regarded 
as a broad distribution and a peak around Q 2 « 0.3 (GeV/c) 2 
not present in the smoother Galster fit. 

In order to get some insight into the consequences of this 
alternative form we have added a term to the form of eq. (1) 
which is able to describe an additional peak with reasonable 
boundary conditions. 



G En {Q 2 ) 



Q 2 



(i + 6Q 2 + c g 4 )5 

dQ 2 

(l + eQ 2 ){l + fQ 2 f 



+ 



(2) 



The rms radius is now given by the sum of a and d, constrained 
to (a + d){2m n ) 2 = a G = 1-73, and we fixed a and d to 
a = 0.37 (GeV/c)" 2 and d = 0.12 (GeV/c)- 2 . The parame- 
ters e and / were kept fixed at 0.5 (GeV/c)~ 2 . Minimising x 2 
yielded b = 0.39 (GeV/c)- 2 and c = 1.68 (GeV/c)" 4 . Here 
we only want to have a parametrisation which reproduces the 
data within the experimental error bars without associating any 
particular physical meaning to the single parameters. In fact, 
as seen in fig. 1, this form reproduces the data well. It is not 
meaningful to go into any detail of an error analysis, instead 
we only show by the example of the dashed-dotted curve that 
with above parametrisation the "peak-region" and the tail to 
higher momentum transfers are essentially described indepen- 
dently from each other. For completeness we just mention that 
the x 2 of the Galster form is by Ax 2 = 4.8 bigger than that of 
the two others. 

As is well known [16], though sometimes questioned (for a 
discussion of this problem see ref. [17]), the Fourier transform 
of the electric and magnetic Sachs form factors Ge{Q 2 ) and 
Gm(Q 2 ) represent the charge and magnetic density distribu- 
tion in the Breit frame, where the energy transfer uj = and 
the three-momentum transfer \q_Breit\ = Q; we denote these 
distributions by p(r), which thus is given by 



p(r) 



4-7T 

1M- 



i: 



'G(QY^l Q 2 dQ. 



(3) 



Refinements to this relation are discussed in detail in ref. [18] 
where it is also pointed out that corrections cannot be defined 
without model assumptions. Since we are interested in the gross 
features of the measured form factors and the spatial distribu- 
tions, we base our further discussion on eq. (3). A more refined 
approach may result in some compression of the resulting dis- 
tributions in r-space, which should not alter their salient fea- 
tures and which are therefore left out of consideration in this 
paper. 

Fig. 2 shows the charge distribution in the neutron, pEn, 
calculated via eq. (3) with above given fits to G_e„. We have 
plotted r 2 pEn{f) which represents the charge in a spherical 
shell at radius r. The charge distribution of the Galster fit shows 
the well known "aperiodic" shape with a positive bump in the 
interior and a negative bump at the outside of the neutron. This 
characteristic feature also results from an ansatz for the form 
factor with the superposition of two appropriate dipole forms, 
to which the Galster parameterisation is a good approximation. 
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Fig. 2. The differential radial charge distribution of the neutron in the 
Breit frame as derived by a Fourier transform. The coding of the lines 
is that of fig. 1. 

The fit with eq. (2), however, which accounts for the bump 
in Geu at Q 2 « 0.3 (GeV/c) 2 , results in an oscillatory be- 
haviour of pEn (r) (see fig. 2). Though the oscillatory behaviour 
depends on the particular fitting form we shall show that it is 
the bump which shifts more charge to the outside than does the 
Galster fit. Since this outer region should be dominated by the 
pion cloud, the corresponding contribution should show up as a 
general feature also in the other form factors, where, however, 
a form factor bump of the same order of magnitude can only be 
expected to be a few-percent contribution. 

With this in mind, we reconsider all four nucleon form fac- 
tors in the following. 



3 The data base 

Table 1 gives an overview of the data which we have taken into 
consideration together with the Q 2 -ranges which they cover. 

For Ge p we have omitted in the final analysis the data by 
Andivahis et al. [22]. In the Q 2 -range of these data, Ge p <C 
G Mp, thus its determination via a Rosenbluth separation is quite 
uncertain. In fact these data are clearly incompatible with the 
new results from polarisation measurements in which not the 
sum of G 2 Ep + t ■ G 2 Mp is measured but the ratio Ge p /Gm p - It 
is straightforward to determine Ge p from this ratio if one takes 
the prevailing Gm p as known from measurement. 

For Gm p we took into account the same data as Kelly [18], 
i. e. the data by Hoehler et al. [29] up to Q 2 = 0.15 (GeV/c) 2 
and those revised and compiled recently by Brash et al. [45]. In 
addition we also used the data by Hanson et al. [23]. 

For Getl we have only taken into account the data from 
polarisation measurements. The measurement in [35,23] give 
only G En , thus the sign of Geu remains undetermined, and the 
errors are so large that the data can essentially be regarded as 
upper limits only; we did not take them into account in the fits. 
Other determinations of Geu were very uncertain due to the 
model dependency of the extraction of G En from the measured 
cross sections, and we did not take them into consideration. 

For Gmu the data by Markowitz et al. [39] and by Bruins 
et al. [40] were omitted in the analysis as was already done 



Measurement 


Q 2 -range 


reference 










p(e,e ) 


0.01 - 


0.05 


f ; . . . nt „i rim 

Simon et al. [19J 




0.04- 


1.75 


Price et al. L20J 




0.39- 


1.95 


Berger et al. |21J 




(1.75 - 


8.83) 


Anaivanis et al. yll\ 


d(e,ep) 


0.27 - 


1.76 


Hanson et al [23] 


p(e,e P) 


0.37- 


0.44 


t** ' 1*1 a i Ai 

Pospischil et al. [24J 




0.38- 


0.50 


Milbrath et al. [25 J 




0.40 




Dietench et al. [26] 




0.49- 


3.47 


Jones et al. [27] 




3.50- 


5.54 


Crayou et al. L28J 










p(e,e) 


0.02- 


0.15 


tt ii ■ i r^m 

Hoehler et al. [29] 




0.16- 


0.86 


Janssens et al. [30] 




0.39- 


1.75 


Berger et al. [21] 




0.67- 


3.00 


cartel et al. [31 J 




1.00- 


3.00 


Walker et al. [32] 




1.50- 


3.75 


Litt et al. [33J 




1.75 - 


7.00 


Andivahis et al. [22J 




2.86- 


31.2 


Sill et al. [34J 


d{e,ep) 


0.27 - 


1.76 


Hanson et al L23J 


Geti 








d(e,e' n)p 


0.15 
0.26 




Herberg et al. [3\ 
bden et al. [1J 




0.30, 0.58 


Seimetz et al. [10] 




0.34 




Ostrick et al. [2] 




0.49- 


1.47 


MaHpv pf al l~1 ? 1^1 




0.76 




filmier et al Till 

vJ ItlZjlCl GL ill • |_J.J_J 




1.00 




Dav pf al TlZLl 
LJay CL dl.|_ l^J 


a {c- ,e njp 


0.21 




r aaSClllCI CL al. [H-J 


He(e , e nj 


0.50 
0.40 
0.67 




7T»n Pt al TQ1 

ZjIIU cl al. \ y J 

Becker et al. [6-8J 
Rohe et al. [5, 8J 


d(e, e ) 


(0.27 - 


1.76) 


Hanson et al [23] 




(1.75- 


4.00) 


Lung et al. L35J 


Gmti 








d(e, e n)p 


0.07- 


0.89 


Kubon et al. [36] 




0.10- 


0.20 


Xu et al. [37] 




0.11 




Anklin et al 1381 




(0.11 - 


0.26) 


Markowitz et al. [39] 




(0.13- 


0.61) 


Bruins et al. [40] 




0.24- 


0.78 


Anklin et al. [43] 


d(e,e'p) 


(0.27 - 


1.76) 


Hanson et al [23] 


d(e,e') 


1.75 - 


4.00 


Lung et al. [35] 




2.50- 


10.0 


Rock et al. [44] 



Table 1. Overview of data taken into consideration (Q 2 in (GeV/c) 2 ). 
The data left out in the final analysis are put into parentheses. The 
reactions are as indicated; <i(e, e') refers to quasi elastic scattering. 
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(and argued) in Kubon et al. [36]. Also, the data by Hanson et 
al. [23] are omitted since they deviate substantially from the 
trend of the more recent data measured with the detection of 
the neutron in coincidence. 1 



4 A phenomenological description of the 
nucleon form factors 

4.1 The phenomenological ansatz 

What we are particularly interested in is the existence of a small 
bump on top of a large "smooth main part". This raises the 
question after what is "smooth main part" and what is "bump". 
After the investigation of several parameterisations we decided 
to keep close to what one is used to in the description of the 
nucleon form factors, namely the dipole form. In fact it is fas- 
cinating that the three form factors Ge p , Gm p , and Gmu are 
describable to quite a precision by the dipole form with the 
one parameter mj given above as it has entered into the text 
books, e. g. [16,47]. Such good description could make one be- 
lieve that there is some physical meaning in the parameter mp. 
In fact there is none. 

Looking at the form factors more closely, however, the pre- 
cision with which the measurements are reproduced by this 
simple parameterisation is limited. While this has been realized 
on the percentage level already in the high-precision measure- 
ment of Ge p at low Q 2 by Simon et al. [19], this became com- 
pletely obvious at high Q 2 by the polarisation measurements 
by Jones et al. [27]. 

Purely phenomenologically, we describe the smooth part 
of the form factors, G S (Q 2 ), by the dipole form. In order to be 
somewhat flexible, however, we took into account the superpo- 
sition of two dipoles: 



G S (Q 2 ) = 



aio 



+ ■ 



«20 



(l + Q 2 / ail ) 2 (l + Q 2 /a 21 ) 



(4) 



To account for a possible bump on top of the smooth form fac- 
tor we now take a parameterisation which is easy to handle 
and the parameters of which give direct insight into the char- 
acteristic features of such a bump, namely its amplitude, po- 
sition and width, ab, Qb, and ab, respectively. A quite natural 
choice would be a Gaussian positioned at Qb- For Qb ^ 0, 
however, such Gaussian contains uneven powers in Q, which 
is not allowed for a function representing a form factor. This 
shortcoming can be healed by superimposing two Gaussians 
as introduced in r-space by Sick for his ansatz for a model- 
independent analysis of nuclear charge distributions [48]. We 
thus parameterise the bump as 



G b {Q 2 ) = e 2V a b 



hi 



Q-Q b \2 



1( Q+Qb \2 

e 2\ a b 1 



(5) 



1 After these investigations were finished a new measurement of 
Gmti was published [46] and came to our attention. Since the Gmu 

values were derived from the reaction s He( e , e'n) they carry a large 
systematic error due to the correction of the binding of the neutron 
in 3 He. They lie even below the values of ref. [43] and support the 
same trend. We repeated our fits with these values taking the error 
into account and found that the results presented in the following are 
not changed. 



In order to keep as close as possible to the accustomed de- 
scription of the form factors, we attribute the full normalisation 
to the dominating smooth part, i. e. Gn(Q 2 = 0) = G S (Q 2 = 
0), i. e. in this ansatz the smooth part accounts for the full 
charge or magnetic moment, respectively. To make sure, that, 
independent of the fitted values for its parameters, Gb does not 
interfere with this normalisation, we multiply it by Q 2 . We thus 
parameterise the nucleon form factors by the ansatz 

G N (Q 2 ) = G s {Q 2 )+a b -Q 2 Gb{Q 2 ), (6) 
where ab is essentially the amplitude of the bump. 

4.2 Fit of the form factors with the 
phenomenological ansatz 

The parameters from the fits of above phenomenological ansatz 
to the form factors Ge p , Gm p , Geu, and Gmu are compiled 
in table 2. The given errors are the standard errors from the fit 
procedure, which also account for correlations, therefore the 
parameters cannot just be varied independently within these 
margins. 

The main purpose of these fits is to allow a coherent view 
on the measured form factors in order to reveal certain common 
features. Therefore, we do not go into details of these fits. How- 
ever, a large number of fits with other analytical forms were 
also tried. They all point to the same feature of a structure with 
a width of « 0.2(GeV/c) 2 . 

The final fitting results given here have been performed 
with fixed normalisation which was guaranteed by setting 020 = 
1 — aio and fitting only aio to the data (a 20 = — »io for Geu)- 
Fits with both a w and a 20 as free parameters did not improve 
the fits by more than Ax 2 = —2; this means that there is no 
hint for normalisation problems in the data. 

Let us first look at the "standard form factors" Ge p , Gm p , 
and Gmu, which we refer to as G s ^ d . For these, one dipole (the 
first) accounts for the overwhelming part of the strength at low 
q, i. e. it carries most of the charge or magnetisation, respec- 
tively. The slope constants of all three G% d essentially agree 
within the errors. In fact, the mean value 0.76 (GeV/c)~ 2 is 
near to the value m 2 D = 0.71 (GeV/c) -2 of the standard dipole 
fit, however the deviation is significant. In fact one cannot ex- 
pect equal parameters since the standard dipole fit deviates sys- 
tematically from the data (c.f. fig. 4) and thus also from our fits 
which reproduce the data within the experimental errors. 

It is interesting to note, that also the slope parameters a 2 i 
of the second dipole form are very similar for all three G s ^ d 
(cf. fit 1 for Gm p in table 2). This term (with negative ampli- 
tude) accounts for the fact that the measurements fall below 
the dipole fit at larger Q 2 . While this became obvious for Ge p 
at larger Q 2 from polarisation measurements, a systematic de- 
viation from the dipole fit was already observed by Simon et 
al. [19] at low Q 2 , though only at the percent level. It may now 
be somewhat surprising that the slope parameters a 2 i are so 
similar for all three G s ^ d . In fact, if one accounts in Gm p also 
for the high Q region (fit 2), this is no longer the case. 

A direct interpretation of the bump structure in terms of the 
parameters ab, Qb, and Ob is obscured by the multiplication 
of the Gaussian with Q 2 . Therefore, for a discussion of this 
structure we refer to its graphical representation in fig. 5. 
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G s 


G b 




form factor 


aio 


an 

(GeV/c) 2 




(221 

(GeV/c) 2 


at 

(GeV/cy 2 


Qb 

(GeV/c) 


(GeV/c) 


Nd.f. 


2 

Xd.f. 


\total 


Ge p 


1.041(40) 


0.765(66) 


-0.041(-) 


6.2(5.0) 


-0.009(7) 


0.07(88) 


0.27(29) 


64 


0.933 


59.71 


GmpI pp 

1 

2 


1.002(7) 
1.003(7) 


0.749(6) 
0.753(2) 


-0.002(-) 
-0.003(-) 


6.0(3.4) 
16.9(6) 


-0.005(1) 
-0.006(2) 


0.35(7) 
0.33(7) 


0.21(3) 
0.23(3) 


60 
75 


0.861 
0.876 


51.66 
65.7 


GEn 


1.04(10.7) 


1.73(-) 


-1.04(-) 


1.54(1.94) 


0.009(6) 


0.29(17) 


0.20(9) 


10 


0.861 


8.61 


GMn/fJ-n 


1.012(6) 


0.770(10) 


-0.012(-) 


6.8(3.0) 


-0.011(1) 


0.33(3) 


0.14(2) 


14 


0.579 


8.11 



Table 2. Parameters from the fit of the phenomenological ansatz eq. (6) to the electric and magnetic nuclear form factors. In the usual way the 
errors on the parameters are given in brackets; if no decimal point is given they refer to the last given digits of the parameter. For Gm p , fit 1 
uses only the data up to 2 GeV/c, thus it is more comparable to the fits to Ge v and Gmu than the fit 2 where all data up to 6 GeV/c are taken 
into account. 



First, fig. 3 shows the overall behaviour of the nucleon form 
factors and the quality of the overall agreement with the fits. 
For Ge p the relatively large (negative) amplitude 020 of the 
second dipole results in a zero in the form factor around 3 GeV/c 
This makes this form factor look very differently from Gmp 
and Gmu, though, in fact, this is only due to the larger am- 
plitude of the dipole with negative sign and not the form of the 
single contributions. At the highest measured momentum trans- 
fers above 3 GeV/c the data for Gm p are not so well described 
by this phenomenological fit. Though, from \ 2 , the overall de- 
scription of the data is excellent, there might be some system- 
atic deviation at high Q 2 pointing to insufficient flexibility of 
the ansatz. 

In order to make the deviation of the measurements (and 
thus also of our fits) from the standard dipole fit more obvious, 
we show in fig. 4 the ratio of the three standard form factors 
(data and fit) to the standard dipole fit (see second factor of 
eq. (1)). 

In fig. 5 we demonstrate the "bump"-contribution to the 
form factors by the subtraction of the fitted smooth part, i. 
e. the two dipoles. In order to emphasise the low-Q 2 region 
where this phenomenon occurs, we have plotted this difference 
as function of \og(Q 2 ). It is obvious from this graph that the 
bump structure around Q 2 = 0.2 — 0.3 (GeV/c) 2 , as discussed 
in section 2 for Gb„, is a common feature of all four form fac- 
tors. It is striking how similar the contribution from the bump 
is in all four form factors. While it is only a small contribution 
to G s ^ d , it dominates G en at low momentum transfer where the 
contribution from the two dipoles, though both separately are 
large, cancel. 

As also discussed in section 2, such bump structure in Q- 
space contributes a certain Apb in r-space, the detailed struc- 
ture of which, however, is not easy to foresee. The net charge 
in Apb under discussion here is zero by construction, therefore 



it must show some oscillation. The wavelength A p of this oscil- 
lation is given by the position of the bump in Q-space, Qb, as 
\ p = 2nh/Qb- The damping of the oscillation is related to the 
relative width of the bump. From Qb ~ 0.45 GeV/c = 2.3 fm -1 
there thus results the wavelength A p = 2.7 fm in agreement with 
the oscillation shown in fig. 2. We conclude that an oscillation 
with such wavelength results as a common feature from all four 
form factors and is not just a peculiarity of G En- 

We demonstrate this by showing in the following subsec- 
tion the Fourier transforms of the four form factors. 



4.3 The Fourier transform of the fits of the form 
factors 

As mentioned above, the Fourier transforms of the Sachs form 
factors can be regarded as the charge and magnetic distribution, 
respectively, in the Breit frame. With this in mind we show 
in this subsection the Fourier transforms of the form factors 
which, for brevity, we all denote by "charge" p. 

Fig. 6 shows p(r) resulting from the two dipole contribu- 
tions, the bump and their sum. For reference, the transform of 
the standard dipole is also shown. In this representation, for all 
three G s ^ d p(r) is very close to that of the standard dipole form 
for r > 0.2 fm. This is true for Gm p also down to the centre 
of the nucleon, where we see deviations from the dipole form 
for Ge p and Gun- It has to be admitted, however, that in the 
latter cases the measurements do not extend to as high Q 2 as in 
the former, thus the distribution in r-space is less well fixed at 
small r. Actually, there is only a tiny fraction of the total charge 
contained in this inner part. The contribution from the bump is 
not visible in this plot. 

The quantity p(r)-r 2 gives directly the weight of the charge 
contained in a spherical shell at distance r, thus the area under 
the curve gives the total charge contained in the respective term. 
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This quantity is shown in fig. 7. Here, the contribution from the 
bump in the form factor is clearly visible as oscillation (net 
charge = 0). Its phase in r-space is such that it puts additional 
strength on the dipole form in the outer region with maxima be- 
tween 1.5 (Gm p ) and 2.0 fm (Ge p , Gmu)- The second dipole 
gives small and tiny contributions in the interior of p(Ge p ) and 
p(Gmu), respectively, and is not visible in p(Gm p )- For Ge u 
the oscillation gives the total p(r) in the outer region centred 
around 1.7 fm, while the inner part is dominated by the dif- 
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ference of the two dipoles. In between there is a cancellation 
between these two terms, a feature already visible in the analy- 
sis in section 2. We come back to this point below. 

In order to emphasise the outer region of the nucleons even 
more, fig. 8 shows \p(r)\ ■ r 2 but now in a logarithmic scale. 
In this plot the sign information gets lost and one has to look 
at fig. 7 to keep track of the sign. The sign of the first lobe in 
the contribution Apb(r) from the bump in the form factor, this 
means abQ 2 Gb{Q), is that of Apb(r = 0). Since Apb(r = 
0) oc a b j Gb(Q)Q 2 dQ, the sign of Ap b (r — 0) and thus the 
sign of the first lobe in Apb(r) is given by the sign of the ampli- 
tude a b of the bump (Gb(Q) > 0). - Here again, we do not want 
to go into details. We will see in the next section, that one ± 
oscillation on top of a smooth distribution can be interpreted as 
dissociation of the nucleon in its counterpart and a pion cloud. 
The further oscillations visible in the logarithmic plot in fig. 8 
are compatible with the data, but certainly must be regarded as 
depending on the special ansatz eq. (6). 



5 A coherent description of the four form 
factors by a physically motivated ansatz 

5.1 The ansatz 

Inspired by the conspicuous graphical representation of the form 
factors, which reveal a bump on top of a smooth trend, we make 
the ansatz of describing the nucleons by the sum of a bare nu- 
cleon plus a polarisation part according to 



p = a p -p 



+ b p ■ (n° + 7T+) 



-p° + n° 



a n ■ n 



+ b p 
+ b n ■ {p u +0 



7T+). 



n a + b n -{+p°-n°+TT-) 



(7) 



(8) 



where we have made use of the normalisation condition ajv + 
6at = 1 for N — p,n. Neutral pions are not taken into con- 
sideration since, to first order, they do not contribute to elastic 
electron scattering. 

The form factors can thus be written as 

G P = G° p + b p ■ {-G°p + G° n + G* + ) = G° p + Gf , (9) 
G n = G° n + b n ■ {+G°p -G° n + G^) = G° n + Gl° l (10) 

where we use the transparent nomenclature of the form factor 
of the polarisation 



G p N oL = b 



N 



(G% - G% + G™"), 



(11) 



where N denotes the neutron (proton) and ir the ir + when 
N is the proton (neutron). 

Furthermore, we think of the bare nucleons in terms of their 
constituent quark content, i. e. p — (uud) and n — (udd). 
We denote the form factor of the distribution of quark q in the 
nucleon N by G qN . We are thus dealing with the ingredients 
G up , G dp , G un , and G dn for which we take the dipole form: 



qN 



G q = 



(i + Q 2 /«D 



qN\ 2 



(12) 
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The pion has intrinsic parity -1 which has to be compen- 
sated by its spatial wavefunction; therefore, to lowest order, it 
should be in an (7=l)-state. Taking as simple ansatz the wave 
function of a harmonic oscillator, the related form factor is 
given by 



= al ■ (1 - l(Q/al) 2 )e-^ Q/ < )2 . 



(13) 



The form factor of the pion cloud should be the convolution of 
this form factor from the wave function with that of the intrin- 
sic distribution of the pion, the size of which is certainly not 
negligible compared to that of the nucleon and thus to the ex- 
tension of the pion cloud. Convolution in r-space results in a 
multiplication in Q-space. Assuming a Gaussian for the intrin- 
sic pion distribution, this results in a multiplication of eq. (13) 
by a Gaussian, thus by a change of the parameter a" in the ex- 
ponential. This does not change the form of G 7r , it would only 
decouple the parameter a\ in the exponential from that in the 
brackets. We will, however, not make use of this additional de- 
gree of freedom. 



For the electric form factors, we have aJ\ 



0.E 





G 


G dp ~ 


G un 


G" r+ ~ 


-G" 



1. For the magnetic form factor, the situation is not that clear 
due to the degrees of freedom of the vector coupling of the 
magnetic moments. Furthermore, it is not clear what magnetic 
moment should be related with the pion cloud since we are not 
dealing with a free pion. 

Strict isospin invariance would imply 



(14) 



We will check, whether the measured form factors can be de- 
scribed under this condition. 



5.2 The electric form factors 

Since G a En (0) = due to the vanishing charge of the neutron, 
from eq. (9) we have 

G Ep (0) = [(1 - &„) • G° Ep (0) + b p ■ G^+ (0)] = 1 . (15) 

We note, that the charge b p ■ 1 of the pion cloud, which goes 
at the expense of the bare proton's charge, contributes to the 
proton electric form factor at Q 2 = 0. Therefore, the peak on top 
of a smooth part of the form factors, which we have revealed in 
section 4, cannot be attributed directly to the pion. In fact, to the 
contrary, according to eq. (13) the contribution from the pion 
cloud should be concentrated around Q 2 = since the pion 
cloud is expected to extend further out than the bare proton. 

According to our model, we evaluate the electric form fac- 
tor of the proton with the ansatz 

G Ep = {G u i + G d l) (16) 
+ b P ■ {-{G u i + Gf) + (G E n + Gf) + Gl + ) 

with G qN and G v parametrised by eqs. (12) and (13), respec- 
tively. 
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The weights for the electric quark form factors are 
given by the quark charges, i. e. a^ p E = +4/3, a^ n E = +2/3 

and a^ p E = —1/3, a dn E = —2/3, and that of the pion is ag + 
■ I This ansatz conserves automatically the normalisation. We 
are thus left with the free parameters of the quark distributions, 
a,i P ,af p ,ai n ,af n , the amplitude b p of the polarisation term 
and the oscillator parameter for the pion, a\ + , i. e. we have 
one free parameter less than in the phenomenological model 
with two dipoles and the bump discussed in section 4. 

In a first fit we take these six parameters as free. The re- 
sulting values are given in table 3 as fit 1 . It is not too surpris- 
ing that the n°-parameters remain completely undetermined. 
The data are described by this ansatz as well as with the phe- 
nomenological ansatz with seven parameters (x 2 = 59.6 here 
compared to 59.7 there). We note that the pion parameters in 
this model are better determined than the bump parameters in 
the phenomenological ansatz. Further, the large values of af 1 
would yield an extremely sharp localisation of the quarks in 
the neutron, most likely an unrealistic scenario. The large er- 
ror on the neutron parameters, however, leave room for ap- 
plying further model restrictions. In fit 2 we subject the quark 
distributions to complete isospin invariance, i. e. we demand 
a"" = a dp and af n = a\ p , there are thus only 4 free param- 
eters left. The proton parameters vary essentially only within 
their errors, the same is true for the pion-cloud parameters. The 
total x 2 increases by the omission of the two parameters by 
only 1.5 which is an insignificant increase. 

The fit 2 is compared in fig. 9 to the data. Here, we also 
show the single contributions of the model. In the logarithmic 
plot (upper panel) their signs get lost, therefore, in order to dis- 
cuss the interplay between the single contributions to Ge p , we 
have plotted in the lower panel these contributions on a lin- 
ear scale. The dominating G° from the bare proton p° has a 
zero around 2.2 GeV/c due to the inference between the posi- 
tive G up and the negative G dp , the latter being suppressed by 
a factor 4, but extending out much further (a dp > a" p ; these 
contributions are not shown separately). This minimum, how- 
ever, is shadowed by the polarisation term G po1 . At Q 2 = 0, 

-b p G° and b p G^ cancel, while G° n itself is zero there (see 
lower panel of fig. 9). With increasing Q, the negative con- 
tribution of -bpGp prevails, G po1 thus becomes negative, until 
it is balanced by the positive contribution from the neutron, 
b p G° n . Around Q w 0.4 GeV/c G po1 has a minimum, around 
Q « 1.0 GeV/c it passes through zero and becomes positive 
at large Q where the contribution 6 p G° prevails. Finally there 
results a zero in Ge p around 3.3 GeV/c due to the interference 
of the positive polarisation (from n°) and the negative lobe of 
Gp from its d-quark contribution. - The data are certainly not 
sufficient to fix these numbers precisely; nevertheless they are 
useful as a guidance for what might go on physically in this 
Q 2 -range. The determination of the minimum in Ge p by ex- 
periment is highly desirable. 

Comparison with the evaluation in terms of the purely phe- 
nomenological model in section 4 reveals that the bump struc- 
ture there has a different meaning than that of the polarisation 
term resulting here from the evaluation in terms of the quark 
model. It is clear, that this is due to the definition of what is 
"bump" and what is "smooth". The interpretation of this model, 
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Fig. 9. Ge p '- Fit 2 (cf. table 3) compared to the measurements. Upper 
panel: Full scale comparison. Lower panel: The contributions to Ge p 
in a linear scale for a discussion of the interplay of the different contri- 
butions to G poL and, finally, the making up of G Ep as the sum of Gp 
and G po1 (note the extended x-scale in the lower panel). - The form 
factors of p° + ty + are normalised to 1, that of n° to 0. 

however, makes clear, that the low-Q 2 side of the bump results 
mainly from the interplay of the form factor from the 7r + with 
the reduction in p°. 

In fig. 10 we show on a logarithmic Q 2 -scale the contribu- 
tion of the pion cloud for fit 2. This scale emphasises the low- 
Q 2 part, where the pion-cloud contribution is concentrated. It 
it obvious that the data are not precise enough to fix the zero 
in G E . A good determination of this zero, however, would be 
a prerequisite to distinguish between the parameter a\ in the 
exponential and in the brackets, thus to see the signature of the 
finite size of the pion. 

We have seen that in the fits to Ge p the parameters of 
the contribution from G" remain practically undetermined. We 
now want to see what we learn about them from the electric 
form factor of the neutron, for which the corresponding expres- 
sion is 

G En = {G U E + G d E n ) (17) 
+ b n ■ {{G up + G dp ) - {GW + Gf) +GI). 
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Fig. 10. The contribution b p G E + from the pion cloud to Ge p (fit 2). 
For comparison the data points show the measurements— [(1 — b p ) ■ 
Ge p + b P ■ G En \. 



As in the case of the proton, this ansatz conserves automati- 
cally the normalisation (here: to 0). Respecting strict isospin 
invariance, i. e. a"" = af p , af n = a\ p ', and G E + — G E we 
can calculate Geu from eq. (17). In the upper panel of fig. 1 1 
we compare this calculation to the measured data, using the 
parameters of fit 2 for Ge p - First, it is remarkable, that and 
how well the polarisation term, directly calculated with the pa- 
rameters from the fit to Ge p , reproduces the bump structure 
of G En at low Q ! However, second, the contribution from the 
bare neutron alone, while being reasonably well positioned in 
Q, overestimates drastically the total measured Geu, 

There are two ways to reduce the amplitude of the super- 
position of two dipoles with equal amplitudes of different sign. 
The first way is to just reduce the common amplitude, here one 
would need a reduction of roughly a factor of 6. The amplitude 
being given by the charge of the one up and the two down 
quarks as +2/3 and -2/3, respectively, such reduction would 
mean to leave the grounds of the present model, namely the 
building up of the nucleons by constituent quarks. Therefore, 
we prefer the second way which requests letting the two pa- 
rameters a"" and af n approach each other. In the fit we let the 
routine search for an appropriate parameter choice by varying 
only a" n , af n , and b p . In fact, the program finds a perfect fit to 
the measured G En with a x 2 per d. f. of 0.81 (fit 2). While fit 
2 reduces the polarisation term by only some 20 %, the dipole 
parameters are increased by a factor of 2.5 and 5, respectively, 
which corresponds to making the distribution in r-space very 
narrow. Furthermore, a™ n /af n =1.2 whereas a^ p /a" p = 2.5, i. 
e. the up- and riown-quark distributions are much more similar 
in n° than in p°. This finding, however, may not be too surpris- 
ing since at the small distances of fractions of a fm the differ- 
ence in the Coulomb interaction in p° and n° might make strict 
isospin symmetry questionable. In other words, the net positive 
charge of the two up constituent quarks will repel them so they 
reside more outside than the quarks with the net zero charge in 
the neutron. One should not mix this up with the opposite be- 
haviour of the current quark distribution as derived from deep 
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Middle panel: a"", af\ and b p fitted to the Gsn-data (see fit 2 to 
G En in table 3). Lower panel: a( n and b„ fitted to the data (see fit 3 
to G En in table 3). 



inelastic scattering. - The result of this fit is shown in the mid- 
dle panel of fig. 11. 

In a last step (fit 3) we examine the significance of the fit 
of a\ n and af n by setting a\ n equal to af from fit 2 of G Ep 
and keeping this fixed. We thus allow only af n and b n to vary. 
While the resulting b n differs by less than one standard devi- 
ation from its value in fit 2, af n just follows a"" in order to 
keep the difference small, which, as said above, is necessary to 
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G 


i 


G 





G n+ 




form factor 


(GeV/c) 2 


of 

(GeV/c) 2 


(GeV/c) 2 


(GeV/c) 2 




ai 

GeV/c 


Nd.f. 


Xd.f. 


Xtotal 


G Ep 

1 

2 
3 


1.000(100) 

1.008(20) 

1.051(20) 


2.03(72) 
2.54(16) 
2.391(14) 


57.(1200.) 

2.54(-) 

2.53(-) 


78.(2500.) 

1.008(-) 

2.22(-) 


0.10(4) 
0.11(1) 
0.118(13) 


0.198(12) 

0.203(6) 

0.204(6) 


64 
66 
66 


0.932 
0.926 
0.928 


59.6 
61.1 
61.2 


Geu 

1 

2 
3 


1.008(-) 
1.008(-) 
1.008(-) 


2.54(-) 
2.54(-) 
2.54(-) 


2.54(-) 

6.2(6.4) 

2.54(-) 


1.008(-) 
5.3(5.1) 
2.22(2) 


O.ll(-) 

0.086(10) 

0.074(5) 


0.203(-) 
0.203(-) 
0.203(-) 


12 
13 


0.807 
0.818 


9.7 
10.6 



Table 3. Parameters from the fits of the model ansatz to the electric nucleon form factors. For the error convention see caption of table 2. The 
different fits observe more or less to isospin symmetry (see text). 





outer 
distr. 


inner 
distr. 


G* 




form factor 


out 

a 


a1 ut 
(Gev/c) 2 


in 

a 


a\ n 

(Gev/c) 2 




a 1 

Gev/c 


Nd.f. 


Xd.f. 


2 

Xtotal 


GMp/l-lp 

1 

2 


0.914(5) 
0.917(6) 


0.818(8) 
0.811(16) 


-0.0049(1) 
-0.0034(14) 


9.578(1.2) 
13.57(6.0) 


0.110(7) 
0.106(8) 


0.213(7) 
0.210(8) 


75 
69 


0.887 
0.901 


66.5 
62.2 


GmtiI i±n 

1 

2 

3 


1.019(14) 

1.363(3.14) 

1.189(1.34) 


0.939(110) 
1.173(700) 
1.060(460) 


-0.112(16) 

-0.511(3.17) 

-0.309(1.38) 


2.37(1.1) 

1.789(2.0) 

1.853(1.8) 


0.219(47) 
0.140(46) 
0.120(40) 


0.152(9) 
0.213(-) 
0.189(19) 


14 
15 
15 


0.629 
0.946 
0.837 


8.8 
14.9 
12.6 



Table 4. Parameters for the fits to the magnetic form factors. In all but fit 3 for Gmu the normalisation is free. Gm p (Qmax - 31.2 (GeV/c) 2 ) : 
Fit 1: Fit to all data, all parameters free. Fit 2: Fit to the data up to 10 (GeV/c) 2 with all parameters free. Gu-a (Qmax = 10 (GeV/c) 2 ): Fit 1: 
All parameters free. Fit 2: af kept fixed at 0.213 GeV/c as determined for Gm p - Fit 3: Normalisation kept fixed by adding a point with value 
1.0000 ±0.0001 at Q 2 =0. 



kept fixed to the values determined now from the fit to the neu- 
tron data. The resulting parameters are shown as fit 3 of Ge p 
in table 3. In fact the change in the parameters and thus also in 
the graphical representation of the form factor are so small that 
we need not go into any detail here. 

Summarising up for the electric form factors, Ge p and Geu 
can be described on the same footing by our constituent-quark- 
pion ansatz. 

5.3 The magnetic form factors 

For Gm p data are measured up to Q 2 = 30 (GeV/c) 2 , there- 
fore in this respect the situation is more favourable here. On the 
other side, the interpretation of the magnetic form factor within 
our model is hampered by the additional degree of freedom of 
the vector coupling of the spins and the magnetic moments: 
While it is clear, that, e. g., the two u-quarks in the proton carry 



keep the Geu small in the high-Q region. As expected from 
the large uncertainties in a™ 1 and af n in fit 2, \ 2 oni y varies by 
0.9 with this quite drastic variation in a"™. The result of this fit 
is shown in the lower panel of fig. 11. 

Data at higher momentum transfers are needed to further 
constrain the low-distance behaviour of the neutron form fac- 
tor. 

We have checked the significance of the bump-structure in 
Geu by fitting the data with only a smooth ansatz consisting 
of two dipoles with equal but opposite amplitudes, which is 
equivalent to the ansatz eq. (1). With this parameterisation we 
get Xtotal = 11-1 (d-f- = 12), i. e. an increase by 1.4 compared 
to fit 2 in table 3. Thus, for a significant determination of the 
bump more precise data at low Q 2 are needed as well as data 
extending to higher Q 2 . 

Finally, we check whether the fit of the proton's electric 
form factor is deteriorated when the neutron parameters are 
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the charge 2-2/3, the resulting magnetic moment depends on 
the coupling of the quark spins. Furthermore, it is not clear 
what magnetic moment one has to attribute to the constituent 
quarks. The same uncertainty holds for the contribution of the 
pion. Though it should predominantly be in a p-state, the re- 
lated magnetic moment is not known, since the pion is highly 
off-mass shell and therefore its mass is not that of the free pion. 
Furthermore, its contribution to the total magnetic moment de- 
pends on the vector coupling. Therefore, in the evaluation of 
the magnetic form factors with the ansatz eqs. (16,17) with 
eqs. (12,13) we have to take also the amplitudes a^ N and 
as free parameters. 

On the other side, one might think that the parameters a\ N , 
which describe the spatial distributions, might be the same for 
Gm and G e such that they can be taken from there. One could, 
however, only profit from this for the sufficiently well deter- 
mined bare proton part, and here only for the dominating term 
from the u-quarks. However, it is not clear whether the dipole 
parameter, determined at relatively low Q 2 , really should hold 
up to the highest Q 2 to which Gm p has been measured. Fur- 
thermore, the magnetic operator does not weight the distribu- 
tion in the same way as does the electric operator. Therefore, 
also the parameters a\ N have to be taken as free. 

Isospin symmetry would suggest that there are only two 
different distributions, that for u- and d-quark in the proton, 
and in the neutron, respectively, the inverted case. In this case, 
including the pion there are only three distributions left and we 
try the ansatz 



Mp'Mp 



G 



M 



Gout , in s~ii' 
+ Cln ■ G 



(18) 



Here, the nomenclature reminds on inner and outer quarks, 
and we omit any discussion about the coupling of their mag- 
netic moments by just giving free amplitudes a™ t,m to their 
respective contributions to the magnetic form factors which, 
again, are parametrised by the dipole form eq. (12) with the 
free parameters a°"*' m . In the same way we allow for a free 
amplitude for the pion cloud. With the a as free parameters, 
normalisation is not guaranteed. - We only mention in passing, 
that we have checked fits with three dipoles; even in the case of 
Gmp, however, up to the highest Q 2 two dipoles are sufficient. 

The parameters from the fits are tabulated in table 4. In the 
fits 1, all 6 parameters of the ansatz were free. Again x 2 is 
comparable to the data evaluation with the phenomenological 
ansatz. We show in the upper panels of figs. 12 and 13 how well 
the data are described. Here, the three terms are also shown 
separately. 

For Gm p we find a surprisingly large value for a™, cor- 
responding to a concentration of the respective distribution in 
r-space near the origin (see subsection 5.4 below), however 
with very small amplitude. For the sake of comparison with 
Gmu, we have repeated the fit with restricting the data to the 
Q 2 -range for which there are data for both magnetic form fac- 
tors (fit 2). We find such large values for ^"(Gjwp) also from 
this restricted data base. About 90% of the (positive) mag- 
netic moment of the proton is carried by the outer distribution 

(ag Mt (G Mp ) ~ 0.91) and 10% by the pion cloud (a5 + (G Mp ) « 
0.11). Note, that the normalisation is violated by some 2%. 
This should be acceptable in view of the quality of the data at 
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Fig. 12. The magnetic form factor of the proton (parameters from fit 
1). Upper panel: Comparison of the measured data with G ou + G %n + 
Gtt + (total) and the three contributions to the fit separately. Lower 
panel: The data points show the measurements— [G ou + G ln ], com- 
pared to the pion cloud a J ■ G n . Note that here the data are shown 
as function of log(Q 2 ) in order to emphasise the low-Q 2 region. 



low Q 2 . The inner distribution contributes only about -0.4% to 
the magnetic moment, while its contribution to the form factor 
becomes comparable at large Q 2 . Further, a° ut (GMp) is 20% 
smaller than a" p (G_E P ), thus the related distribution in r-space 
extends further out for the magnetism than for the charge. 

In fit 1 of Gmu, the dominant contribution to the magnetic 
moment again comes from the outer distribution (note that by 
referring to Gmu/ the signs are inverted). The sign of the 
inner distribution again is negative, in this case, however, its 
contribution to the magnetic moment is about 1 1 % and thus not 
negligible. The pion cloud contributes a factor of two more to 
/i n than to u p in this fit. It has to be admitted, however, that in 
this fit the normalisation is off by about 10%. In fact, the data 
do not extend sufficiently far down in Q 2 to let the normali- 
sation free in the fit, and the pion cloud is particularly sensi- 
tive to the data at low Q 2 . Fit 2 shows the result of a fit with 
af{G M n) fixed to af(G M p) = 0.213 GeV/c. The fit now 
obeys normalisation to within a percent with the amplitudes 
dQ Ut and a " having very large (correlated!) errors. With fit 3 
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Fig. 13. Same as fig. 12, now for Gmti\ fit 3 in table 4. 

we went one step further by adding an additional data point at 
Q 2 = in order to fix the normalisation, while at the same time 
letting the parameter a\ (Gmu) free. There is some redistri- 
bution between inner and outer distribution, but all changes 
of the parameters are within the errors. Thus, there is no prob- 
lem with the normalisation of the data. - We only mention in 
passing that taking into account also the data by Markowitz et 
al. [39] and by Bruins et al. [40] yield df (G Mn ) between 0.09 
and 0.12 and of (G Mn ) between 0.186 and 0.189 GeV/c. 

It is beyond the scope of this analysis to try an explanation 
of these findings. 
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Fig. 14. r 2 ■ p(r) in the Breit frame, calculated with parameters given 
in tables 3 and 4 (Ge p ' Fit 3, Gm p '- Fit 2, Gun'- Fit 3). 



5.4 The distributions in r-space 

In fig. 14 we again show the distributions r 2 ■ p(r) in the Breit 
frame for the three standard form factors, now for the model 
evaluation calculated with parameters given in tables 3 and 4. 
In the proton, the contribution from the inner distribution is 
practically invisible. This shows that, to the degree of precision 
visible in this plot, the proton form factors are describable by 
one dipole plus the contribution from the ir + , which builds a 
shoulder on the distribution extending out beyond 2 fm. The 
magnetic distribution in the neutron has an appreciable con- 
tribution from the inner distribution. Note that by evaluating 



Gmu/ Mn a U signs are inverted such that, e. g., the contribution 
from the tt~ comes in with a positive sign. It is worth to men- 
tion, however, that it is the fit which yields the positive sign 
for the contribution parametrised as form factor of a lp wave 
function. 

To emphasise again the smaller contributions and thus in 
particular the outer region, fig. 15 shows r 2 ■ p(r) in logarith- 
mic scale. By construction, the distinct structure at the edge 
of the distribution now consists of only one bump, which, ac- 
cording to the model, is due to the pion cloud. This evaluation 
shows that the oscillations in the phenomenological analysis in 
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Fig. 16. The contribution of the polarisation term to r ■ p(r) for the 
proton (in the Breit frame). 



bp ■ 7T+ yields just the two lobes which were also seen in the 
phenomenological analysis and which are emphasised in the 
logarithmic representation in fig. 8. The small shift in the zero 
compared to fig. 8 is due to the difference in what is regarded 
as the smooth part of the form factor. In fig. 15 the negative 
lobe of the "bump" in fig. 8 in is not visible since — b p ■ p° is 
absorbed in the contribution from p° as a whole. 

The charge distribution of the neutron, see fig. 17, is dom- 
inated by the smooth polarisation oscillation, i. e. by the posi- 
tive lobe b n ■ p° and the negative lobe from the tt~ . These two 
contributions add up to the same form as the polarisation in 
the proton, but with opposite sign. Superimposed is now the 
charge distribution from the neutron, (1 — b n ) ■ n°, which mod- 
ifies the smooth oscillation in a characteristic way. In partic- 
ular it reduces the positive lobe from b n ■ p° around 0.5 fm, 
possibly leading to a region with zero net charge. The details 
depend on differences in the ansatz, the present data do not con- 
tain sufficient information to discriminate between the different 
solutions. It is, however, gratifying to note that this feature is 
present throughout the different approaches in this paper: It is 
also visible in figs. 2 and, particularly clearly, in the lowest 
panel of fig. 7. 



Fig. 15. Same as fig. 14, now in logarithmic scale. 



section 4 are not significantly determined by the data, they re- 
sult from the particular phenomenological ansatz used there for 
the separation between a "smooth" and a "bump" contribution 
to the form factor. The shoulders in all three standard form fac- 
tors, however, emerge in both evaluations, and we judge them 
as being an unambiguous result from the data. 

In fig. 16 we show the polarisation contributions to the elec- 
tric form factor of the proton in an enlarged scale for a closer 
comparison with the situation in the neutron which is shown 
in fig. 17. The (tiny) neutron contribution to the polarisation 



part of the proton, b p 



is situated in the inner region. The 



superposition of — b p ■ p°, i. e. from the reduction of p°, and 



6 Conclusion 

It is found as a common feature of all four nucleon form fac- 
tors that they exhibit a very similar structure at small momen- 
tum transfer, which is related with some structure in r-space at 
large r around 2 fm. Such finding asks for a common explana- 
tion. We propose to interpret this as resulting from a pion cloud 
around the bare nucleon. This is actually an old idea account- 
ing for the chiral symmetry in quark bag models of the nucleon 
and used since many years [49,50]. The phenomenologically 
successful "cloudy bag model" (see e.g. [51] and references 
therein) was recently used to describe the form factors of the 
nucleon [52]. However, this description was still based on the 
old data base and did not look for the effect of the pion cloud 
in the form factor at low Q 2 . 
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Fig. 17. r 2 ■ p(r) in the neutron (in the Breit frame). The contribu- 
tions from the bare neutron and proton and from the pion are shown 
separately. Upper panel: Fit 2, lower panel: Fit 3 in table 3. 



Also the other ingredient of our polarisation model has a 
deeper theoretical basis. The division between the bare pro- 
ton and neutron contributions into separate contributions from 
up and down quarks is not only suggested by the naive con- 
stituent quark model but is well justified through calculations 
in quenched lattice QCD [53,54]. 

A fit of the data with the ansatz of a p-wave for the pion 
yields probabilities for the dissociation of the nucleons into 
their counter part and a charged pion of 7 to 20%, details being 
dependent on the peculiarities of the ansatz. One could com- 
pare these probabilities with similar results from high energy 
experiments. But without a more thorough theoretical discus- 
sion such a comparison is not very meaningful since our polar- 
isation model is rather crude. It is evident that next-to-leading 
order contributions as the two pion continuum or the pion-Delta 
component in the nucleon wave function have to be considered. 
In a dispersion-theoretical analysis such contributions can be 
included naturally and, in fact, have been in the past [55]. Such 
analyses will be pursued again in the frame work of chiral dy- 
namics particularly suited to take higher order effects to the 
pion cloud into account [56]. 

From the experimental point of view our exciting result 
shows that further studies at momentum transfers squared Q 2 < 



1 (GeV/c) 2 down to the lowest reachable values are much needed 
with increased precision. 

We have parametrised the smooth part of the form fac- 
tors by the superposition of dipoles, which lend themselves 
to an interpretation in terms of the distribution of constituent 
quarks. Data at high momentum transfers are needed to check 
this model assumption and constrain the distribution. 
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